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Abstract 

For a self-repelling polymer chain consisting of n segments we calculate the persistence length 

^ ; 

^ . L{j,n), defined as the projection of the end-to-end vector on the direction of the j^^ segment. 

, This quantity shows some pronounced variation along the chain. Using the renormalization group 

and e-expansion we establish the scaling form and calculate the scaling function to order e^. 

■ Asymptotically the simple result L (j, n) ~ const (j (n — j) /n) ^ emerges for dimension d = 3. 

o , 

C^, Also outside the excluded volume limit L(j,n) is found to behave very similar to the swelling 

' factor of a chain of length j (n — j) jn. We carry through simulations which are found to be in 

■ good accord with our analytical results. For d = 2 both our and previous simulations as well as 

(N : 

. theoretical arguments suggest the existence of logarithmic anomalies. 
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I. INTRODUCTION 



In solution a long flexible polymer chain takes a random, coil-like configuration. Many 
properties of these coils adequately are described by a simple model, where the chain is 
taken as a linear sequence of n structureless segments Sj = rj — rj_i, with vectors 
{j = 0, . . . , n), giving the positions of the endpoints of the segments in (i-dimensional space. 
In such models the chemical microstructure of the polymer is absorbed into a few parameters 
like the mean squared segment size i"^ ~ (s^) or the excluded volume uq, which measures 
the interaction among the segments. For uq > long chains are swollen compared to a 
noninteracting random fiight chain, and it is this 'excluded volume region', which will be 
considered here. We will discuss the infiuence of the excluded volume on the persistence 
length, which measures the range over which the chain configuration on average remembers 
the direction of a specific segment. It is defined ^ as the projection of the end-to-end 
vector r„ — Tq on segment vector s^-: 

„, = . (1.1) 

Here the pointed brackets denote the thermodynamic average. We note that often the 
persistence length is identified with L{l,n), i.e. with the projection on the first segment. 
We here use the more general definition (1.1), since the position along the chain of the 
distinguished segment Sj will turn out to be an essential variable. Indeed, L{j,n), varies 
considerably along the chain. Close to the chain ends {j/n — or 1, respectively) it is 
a nonuniversal, microscopic quantity, whereas in the critical limit {n — > oo, < j /n < 1, 
fixed) it becomes a universal macroscopic observable of the excluded volume coil. 

A microscopic, nonuniversal persistence length is an important parameter of the model of 
a 'worm like' chain l|, where bond-angle constraints correlate the directions of subsequent 
segments, but no excluded volume interactions among segments spaced a large distance 
along the chain exist. In that model the persistence length is a measure of the local stiffness 
of the chain and asymptotically becomes independent of the chain length n. The concept 
of a persistence length also plays a prominent role in theories of polyelectrolyte solutions. 
Here the persistence length refers to an underlying worm-like chain model, and various 
definitions including also L{j,n), Eq. (1.1), are used. (See Refs. for recent work 

and a compilation of literature relevant to this topic.) All this work treats the persistence 



length as a more local quantity, independent of chain length. Excluded volume effects are 
neglected. 

If we take the excluded volume into account, we easily see that the persistence length, as 
defined in Eq. (1.1), cannot be independent of chain length. Rather it must show power law 
behavior as function of n. This follows from a sum rule, relating n) to the mean-squared 
end-to-end distance 

Rl{n) = {{v^-v,f) . (1.2) 
Taking (s^) = const, which should be a very good approximation, we immediately find 

n 

Y,L{j,n)^Rl{n)^n"' , (1.3) 
i=i 

implying 

L(j,n)~n2^-i . (1.4) 

This result should be valid in the excluded volume limit of long self-repelling chains. Since 
u > 1/2 {u ^ 0.588 for d = 3), this shows that L{j,n) is a critical quantity, diverging with 
increasing chain length. It thus is somewhat surprising that it has not found much attention 
among workers concerned with the excluded volume problem. We are aware of only three 
papers , where L{j = l,n) is calculated for self avoiding lattice walks, using exact 

enumeration or Monte Carlo methods. Most results given there are for two-dimensional 
la..., whe. the results st,,est tKat .,1, „) d,ve,es loga.tK.leaUy fl Q o, wltK a ve., 
small power 5] of n. On the cubic lattice L(l, n) seems to tend to a constant [7]. In view of 
the sum rule (1.3) these results immediately imply that L{j,n) must be strongly dependent 
on j, as stressed above. 

In the present work we use renormalized perturbation theory to calculate L{j,n) for a 
chain with short range self-repulsion. We concentrate on the dependence on segment index 
j, where in three dimensions we will find a surprisingly simple asymptotic behavior. In two 
dimensions, however, the behavior of L{j,n) is more complicated, and our results support 
the existence of logarithmic anomalies. For reasons explained later we in our calculation go 
to second order in the renormalized coupling constant (two loop). We furthermore compare 
our results to simulation data for chains up to length n = 2000 on a cubic lattice. Simulation 
results for a square lattice are also given. Our analysis is an extension of our previous work 
on the correlation function (sj^ • Sj^) of individual segment directions. 



This paper is organized as follows. In Sect. |n]we define the model and introduce the 
general structure of perturbation theory. In Sect. IIIll we evaluate L{j,n) to two loop order. 
The one loop result is analyzed with the help of a crossover formalism suggested by a blob 
picture. Such an approach is known [9j to yield good results for a large variety of observables 
for dilute or semidilute polymer solutions. The two loop calculation, evaluated in strict e- 
expansion {e = A — d), serves to support the assumptions implicit in the crossover model. In 
Sect. II VI we compare our results to our simulations in three dimensions and also discuss the 
logarithmic anomalies showing up in (i = 2. Sect. IVl summarizes our results. 



II. THE MODEL 

Our model has been presented in detail in Refs. 0, and we here briefly recall the 
essential features. 

As mentioned in the introduction, we describe the chain configuration by the set of vectors 
Tj {j = 0, . . . ,n). The energy is written as 

H 



V0 + V2 , (2.1) 



where 



i=i ° 

incorporates the connectivity of the chain, and 

l[[l-{Anilf'pj'{r,-r,>)] (2.3) 



-V2 
e ^ 



j<j' 

represents the excluded volume interaction among the segments. Here the excluded volume 
is written as uq = (Anil)'^^'^ Pe, introducing the dimensionless excluded volume parameter Pe. 
The microscopic length scale io governs the segment size. In relation to more microscopic 
modells it incorporates all microstructure effects usually adressed as intrinsic stiffnes of the 
chain. For /3e = one finds 

{s^)o = 2dil . (2.4) 

In Eq. (2.3) the prime at the product indicates that in multiplying out we omit all terms in 
which some vector occurs more than once. 
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To calculate L (j, n) we define a generating functional 

Z(q,h) = ^^""^^ ' j e-Vo-V2gh.s,+»q(r„-ro) ^ ^2.5) 

I?[r]=ff^^!^ , (2.6) 

where f2 denotes the volume of the system. The normalizing factors are chosen such that 
the partition function 

Z(n)=Z(0,0) (2.7) 

reduces to Z{n) = 1 for a noninteracting system {(3e = 0). From 2 (q, h) we can calculate 
the (not normalized) persistence length as 

L (j, n) = (s, ■ (r„ - ro)) = -f- {-t V, ■ V,) Z (q, h) . (2.8) 

Z [n) q = = h 

The generating functional (2.5) is evaluated by expanding the product (2.3) in powers of 
(3e- The individual contributions can be represented by diagrams in which the polymer is 
drawn as a straight line and the excluded volume interaction is represented by broken lines 
('vertices'), connecting pairs of special points which correspond to the interacting segments 
(ja, jp) (see, e.g. Fig.^J. Parts of the polymer line connecting two consecutive special points 
will be addressed as propagator lines, labeled by some 'momentum' variable k. Momentum 
conservation holds for each vertex. The chain ends 0, n represent special points closing 
propagator lines of momentum q or — q, respectively. 

In evaluating a diagram each broken line stands for a factor — {ATriQ)^^"^ j3e, and a propa- 
gator line of momentum k connecting special points ja < jp stands for exp (— k^^Q (j/j — jo-)). 
Internal momenta are integrated over all space: 



and the segment labels of the special points are summed from 1 to n — 1, respecting their 
ordering along the chain. 

With these rules we can evaluate Z (q, 0). The derivative V/i|o Z (q, h) involves a further 
vertex, drawn as a stroke through segment j in the polymer line. It arises from the differ- 
entiation V/i|o and yields a factor 2i^1\^, with k being the internal momentum flowing into 
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• — ^ 

a) 



b) 

FIG. 1: Classes of diagrams discussed in the text. The grey blobs stand for any number of 
interactions. 



Z{n) = 1+ + 



a) 



^ 9 + q + ki^q + ki + k-i N if + ^ ^ , 

31 j 32 33 34 

b) 

FIG. 2: a) Diagrammatic expansion of Z (n) b) Diagrams contributing to Z (n) Z,'^''^'^^ 

that special point. The segment index j is not summed over. An explicit example for the 
application of these rules is given below. 

In evaluating L {j, n) (Eq. 2.8), the contribution of each diagram of the general structure of 
Fig. ^ vanishes. Here Sj is part of a polymer loop closed by the explicitly drawn interaction. 
The loop does not interact with the remainder of the chain and therefore the direction of 
Sj is not correlated with r„ — fq. Diagrams of the structure shown in Fig. ^ yield a very 
simple contribution. Due to momentum conservation the momentum flowing into the stroke 
is q, and under the operator — iVg|o the total contribution of all such diagrams yields 
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2dil Z (j) Z {n — j). We thus find a 'reducible' contribution 

L(-^)(,,n) = 2rf£^^i^l|/!^ . (2.9) 

Z[n) 

Up to order I3l the diagrams for Z{n) are shown in Fig. The remaining 'irreducible' 

contributions L*^''''') (j)"^) are at least of order /?g, as shown in Fig. Et". To exemplify the 

application of the evaluation rules we write down the contribution to L^"'''^ (j, n) of the 
explicitly labeled diagram: 

ii=l 3<32<ji<h<n ^2 

g-q'^giig-(q+ki)'^g(j-ii) _^ j^^) e-(q+ki)'£g02-j) 

g-(q+ki+k2)^£g03-i2)g-(q+k2)^^g04-i3)e-q^^g(n-i4) _ 2rf£2/32Di(j) (2.10) 
ii=l i<i2<j3<j4<n [(j4 - J2) (j3 - Jl) - (j3 - J2) J 

The second diagram of Fig. |2|d yields 2dilP'^D2{j), where 

0<ii<j2<j j<j3<j4<n [{ji - 32) {j3 - Jl) - (j3 - J2) \ 

whereas the last diagram results from the first one by refiection of the chain {j ^ n — j). 
We thus find 

L(-) (j, n) = 2MlPl [D, (j) + D, (j) + D,{n- j)] + O ) . (2.13) 

Now it is clear that only the irreducible diagrams yield new contributions specific for the 
persistence length. It is for this reason that we evaluated L{j,n) including order 



III. CALCULATION OF L {j, n) 

A. First order unrenormalized perturbation theory 

In order (3^ only the first diagram in Fig. contributes: 



0<jl<j2<n 



k 



l-/5e E 02-Jl)-'/VO(/3.^) . (3.1) 
0<ji<j2<n 



Substituting this result into Eq. (2.9) for U^^ (j, n) we find 

L {j, n) = 2del [1 + (3,R, (j) + O {(3!)] , (3.2) 

j-l n-l 

= E E -Ji)-'^' . (3.3) 

jl=l j2=j + l 

For e = A — d > the summations in Ri can be approximated by integrals, which implies 
that we take the limit of a continuous chain. This yields 

= 1 + ^^^^^'7(^ [^'^' + (1 - ^y^' - + ^ (^3 ' (3.4) 

where we introduced the notation 

J=- • (3.5) 
n 

The result (3.4) is correct for n ^ 1, up to terms of relative order 1/n neglected in the 
continuous chain limit. 

B. Renormalization 

Renormalization exploits the fact that for n ^ 1 the microstructure becomes unimportant 
and physical observables are invariant under a change of io, compensated by an appropriate 
change of (3^ and n. The theory has often been explained in the literature, and we here use 



the formulation of Refs. [8|, | 

We introduce a renormalized length scale 



'i , (3.6) 



and we define the renormalized coupling u and the renormalized chain length n^j by the 
formal relations 

= \'uZ^{u) , (3.7) 
n = X~'^nRZn {u) . (3.8) 

The ratio J = j/n is invariant under renormalization. The parameter A obeys the inequality 
< A < 1, but otherwise is arbitrary. The renormalization factors Zu and Zn are chosen to 
absorb the poles in e, which show up in expressions like Eq. (3.4). In the continuous chain 
limit implicit in the evaluation of segment summations as integrals, these poles carry the 
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information on the microstructure. Within the scheme of minimal subtraction the renor- 
mahzation factors up to the order needed here are found as 

Zu{n) = i (^l + ^u + 0{u^)^ , (3.9) 

Zniu) = l---(^-Pju' + 0{u') . (3.10) 



e V2e^ 

The renormahzation factors have been calculated to higher orders in u, and very accurate 
expressions for the logarithmic derivatives dlnu/dln X, din Zn/dln X have been derived, 
known as renormahzation group flow equations (see e.g. Ref. |10|])- For X —* 0, which 
corresponds to the excluded volume limit of long self-repelling chains, the renormalized 
coupling tends to a fixed point u*. With our convention one finds 

M* = 0.364, (rf = 3) . (3.11) 

We also will need the e-expansion 

Integrating the flow equations one derives the renormahzation group mapping from bare to 
renormalized parameters. In terms of a normalized coupling 

7/ 

/ = - (3-13) 

u* 

it reads 

iR = f\l-fr'/^Hu{f)se , (3.14) 
UR = f-'\l-fr^'-^§^s^n . (3.15) 

[J) 

Here si = ioSe{Pe), Sn = s„(/3e) are integration constants, which contain the dependence on 
the bare parameters io, (3^. In practice they are to be taken as microscopic fit parameters. 
uj and u are critical exponents, which in three dimensions take the values 

0.588; a; « 0.80 . (3.16) 

We below also will need the e-expansion of u: 
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Finally, the functions Hu (/) and H (/) within the accuracy of the mapping can be param- 
eterized as 



H^U) = (1 + 0.824/)°-'^ , 
H{f) = 1 - 0.005/ -0.028/2 + 0.022/3 



(3.18) 
(3.19) 



So far we just recalled some general results of renormalization group theory. We now 
turn to the renormalization of the persistence length. The sum rule (1.3) indicates that 
nL {j,n) ~ ((r„ — fq)^). Since Rl (n) = ((r„ — fq)^) is an observable, invariant under the 
renormalization group, this suggests to define the renormalized persistence length as 



n 



Lr (j, ur) = — L (j, n) = X Zn (u) L (j, n) . 
riR 



Indeed, using Eqs. (3.4), (3.6) - (3.10), we find 



Lr (j, ur) 
2dil 

For e — > this yields 

Lr (J, Ur) ^ 
2dil 



— h urin — 

e ^ e(2-e) 



J^/' + (l-j) 



(3.20) 



+ O {u^) (3.21) 



1 + 



u 



l + lnnii + lnJ+ln(l-j) + 0(e) + C (u^) 



(3.22) 



The pole in e is cancelled, as expected for a properly renormalized quantity. 



C. Crossover analysis of the first order result 

Though by construction the leading long chain behavior of observables like or nL (j, n), 
if evaluated to all orders of perturbation theory, is invariant under renormalization, low order 
perturbative approximations generally depend on our choice of the renormalized length scale 
Ir. An exception from this rule is provided by so called 'universal ratios', i.e. dimensionless 
ratios of observables, constructed such that all explicit renormalization factors drop out. 
Evaluated at the fixed point u*, such ratios have a unique e-expansion, independent of any 
conventions of the renormalization scheme. In the next subsection we will construct such 
a ratio. Here we are concerned with the direct evaluation of the result (3.21), and we thus 
have to face the problem of the proper choice of Ir. 

Previous work on many different experimental observables shows that we can con- 
struct a good approximation by evaluating first order results like Eq. (3.21) directly in three 
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dimensions (e = 1), with taken to be of the order of the smallest length relevant to the 
problem considered. In other words, we use the renormalization group to map the physical 
chain on a chain of effective segments ('blobs') of size chosen such that the quantity con- 
sidered does not resolve the internal structure of a blob. For the present problem the blob 
should be identified with the smaller one of the two subchains (0,j), {j,n). To construct a 
smooth crossover among the limits j ^ and j ^ n we define a chain length variable 

. ^ J(r^-j) ^ (3.23) 

n 

and we implicitly fix iji = i^/X by choosing 

nn = X'Z-'h = no , (3.24) 

n 

where no is a constant of order 1. To stay consistent with our previous work 9] we take 
hq = 0.53, a value determined from an analysis of the interpenetration ratio, which is 
proportional to the second virial coefficient of the osmotic pressure divided by the coil 
volume. 

Taking e = 1 we find from Eqs. (3.20), (3.21), (3.23), (3.24) 

L{j,n) = ^LRij.rLR) 
n 



6^4 



n 

Eqs. (3.14), (3.15) yield 



no = r'\^-f\^^Mrv'nn , (3.26) 



f\^^-^-')H{f)sisn . (3.27) 



As a result we find the crossover form of L {j, n) , evaluated to first order renormalized 
perturbation theory: 



L{j,n) = Qs^si\l-f\-^--')Hif) 



1 + 2u*fnl^' (1 - J)-'/' + r'^' - r'^' (1 - J) 



1/2 / n _ , --1/2 _ --1/2 _ -.-1/2 _ 1 -1/2 



n 







(3.28) 



2 

This result deserves a detailed discussion. 

Eq. (3.28) involves the parameter se, which has dimensions of a length, only in the 
combination 

P = Sns] . (3.29) 
11 



This is a general feature of the renormahzed theory. At the 0-point, which in the present 
theory corresponds to a strictly noninteracting chain (/3e = = /), d reduces to £o- (For 
T > 6 it is weakly temperature dependent: i = £o {1 + O (Pe)), Pe ^ T — 6.) Evaluating 
Eq. (3.28) at the 0-point / = we thus find the expected result: 

L(j,n) = 6£^ = (s2), T = e . (3.30) 

The persistence length is of microscopic size, independent of n and j. Note that keeping 
residual three body interactions and treating the 6-point as a tricritical point we in three 
dimensions expect to find logarithmic corrections similar to those found for the end-to-end 
distance ^]: L (j, n) ^ 6£q (1 + const/ Inn). 

The excluded volume limit /5e > 0, n — > oo, is reached for / ^ 1. Using Eq. (3.26) to 
eliminate |1 — / | we from Eq. (3.28) find 



1 + 2u*n'J' ( (1 - J)-^/^ + J-V^ - j-^' (1 - jY 



-1/2 -L -1/2 

nn 

2 ° 



(3.31) 



where, as usual, the star indicates the excluded volume limit. The microscopic length 
parameter B is defined as 

B = s,sy,-'Hl''^m''a) ■ (3.32) 

It is the only microscopic parameter showing up in the excluded volume limit. 

The result (3.31) shows the expected overall scaling L* ~ n"^^'^. It furthermore predicts 
some pronounced dependence of L* (j, n) on J. For J <C 1 it yields 

t (j, n) = 6B'no {nj^'' [l + u* (2nf -l + O {3"^)) 

~ f -^'^ . (3.33) 

Thus the persistence length is microscopic close to the chain ends and rapidly increases for 
J approaching the center of the chain. We also note that L* {j, n) for j <^n (or n — j <^n) 
is independent of the chain length. This feature is quite plausible. It implies that the blob 
(0,j) essentially is infiuenced only by its neighboring blobs along the chain. Starting from 
this hypothesis we easily can derive the basic structure of Eq. (3.31) from a simple scaling 
argument. Assuming that L{j,n) is renormalizable, we from the sum rule (1.3) find the 
scaling law 

t{j,n)=n"'-'C{j) . (3.34) 
12 



Then assuming that the hmits n oo, j fixed, and n oo, {n — j) fixed, exist, we 
immediately find 

C{j) = {j{l-j)r-'C,{j) , (3.35) 

where £i (0) = Ci (1) takes some finite value. 

The variation of L* {j,n) (Eq. (3.31)) closely resembles the behavior of the end-to-end 
swelling factor a"^ [n) of a chain of length n = nj{l — j). The latter is defined as 

«l(n) = ^ , (3.36) 

where the denominator is the mean squared end-to-end distance of a noninteracting refer- 
ence chain, reducing to the physical chain only at the B-point. (Recall the remark below 
Eq. (3.29).) The close similarity among L (j, n) and a'^ (n) holds also outside the excluded 
volume limit. Within our renormalization scheme we to first order renormalized perturbation 
theory find 



1/2 / 2 _i/2 



1 + u*fnQ' [--n 



(3.37) 



From Eqs. (3.28), (3.29), (3.37) we may construct the ratio 
L (j, n 



) 1 + 2«*/nf ((1 - J)-^/^ +J-'/'-r'/' (1 -J)-'^' 



6^^«|(n) l + «*/nf (l-n-^/^ 



(3.38) 



Here the prefactor |1 — f\^^~^', which contains the dominant variation, has dropped out, 
leaving only some weak dependence on J and the coupling strength /. For ] = \ the ratio 
(3.38) varies from 1 for / = to about 1.3 for / = 1. We thus find L {j,n) ~ a"^ (n) quite 
generally. Indeed, if evaluated at the fixed point / = 1, 

_ L{j,n) fiL {j,n) 
' - QPal in) - Rl in) ^'-^'^ 

is an universal ratio, which will be calculated to order in the next subsection. We there 
also will compare the crossover of L [j, n) from 6-conditions (/ = 0) to the excluded volume 
limit (/ = 1), as resulting from the present analysis, to the result of second order e-expansion 
(see Fig. IH). 
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D. Expansion to order 



In the analysis of the previous subsection we imphcitly have built in the power law 
behavior (3.35) 

by choosing the renormalized length scale in to be of the order of the end-to-end distance 
of the smaller subchain. This amounts to incorporating the blob picture underlying the 
scaling approach, and is adequate only if the limit n — oo, j fixed, is finite. The validity 
of that approach is not evident a priory. Other observables related to the internal structure 
of an excluded volume chain, like the distribution of distances among segments ji and j2, 
are known to show specific end effects, showing different power laws close to a chain end as 
compared to the center of the chain. It therefore is conceivable that also L{j,n) shows new 
singularities for j 0. In Sect. IIV Bl we will show this to be the case in two dimensions. 
Here we check the behavior of L*{j,n) within the framework of strict e-expansion. The 
logarithmic terms showing up in this expansion must sum up to the expected power law. 

It is easily checked that the first order result (3.22) obeys this criterion. Using the e- 
expansions of u* (Eq. (3.12)) and u (Eq. (3.17)) we find 

^^^^^ = l + |(l + lnn« + lnj+ln(l-j)) + 0(e^) 
= {l + l){nRj{l-j)f' + 0{e') 

= [l + "-){nnj{l-j)r~' + 0{e') . (3.40) 

However, as pointed out at the end of Sect. IHl the irreducible contributions specific to 
L (j, n) occur first in second order. It thus is appropriate to calculate Lji (j, n) to order e^. 

Unrenormalized expressions for the irreducible diagrams have been given in Sect. ITU 
Eqs. (2.10) - (2.11). For the reducible contribution we find 

= 1 + PeRl (j) + Pi [Rl (j) + i?2 (j) +R2{n- j) 

+ R,(j) + R,(n-j)-R,ij)) + 0{p!) , (3.41) 

where Ri (j) has been defined in Eq. (3.3), and 

i-i 

^2 (j) = E E (-^3 - J2r'^' [(J4 - Jir"^' - (J4 - J3 + J2 - Jl)-'^'] ,(3.42) 

ii=l j<j2<j3<ji<n 
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i-1 



^3 (j) 



R4 U) 



ii=l j<j2<j3<j4<n 



-d/2 



((is - jl) (j4 - i2) - (is - j2f) '^''^ 



Y Y [ (-^3 - h) '^''^ iu - is + i2 - ii) 



-<i/2 



0<Ji<j2<j j<j3<j4<n 



+ ((is - ii) (i4 - i2) - (is - i2) ) 



2\ -d/2 



(3.43) 



(3.44) 



Here the diagrammatic contributions have been combined such that for e > all summations 
can be evaluated as integrals. The evaluation of expressions (3.42) - (3.44), (2.11), (2.12) is 
straightforward, but lengthy. The resulting unrenormalized expansion reads 

L (i, n) 



2dil 



= 1 + f3,n^/'a, (J) + Pln^a2 (j) + O (/3|) , 



with coefficients given in e-expansion as 



(3.45) 



ai(j) = - + l + lnj+ln(l-j) 



«2 (j) 



+1 (^l + In j+ In (1 + l In' J+ ^ (1 - j)^ + O (e') , (3.46) 



6 1 A3 
-^--h7 + 61nj + 61n(l-j) 



e V 2 
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-| + y - f InJ- y ln(l - J) - 31n' J- 31n' (1 - j) 
+ lnjln(l-j)+^(j)+^(l-j) + C»(e) , 



(3.47) 



with 



^(J) = ^ In J In ( 1 - J) - ^4^- In 

4 2V1-J 



In (7 - 9J+ 3/) In' ^"^^"^^ - In' V^- j 



y dt{9- 6t) In' 



y/4:-3t + Vi 

2^/r^ 



(3.48) 



The function (j) contains no leading logarithmic terms, but vanishes for j ^ and stays 
finite for j— > 1. We now use Eqs. (3.6) - (3.10), (3.20) to find the renormalized expression 



2M\ 2 



l + lnfiR + l fl + lnnR + ^ln'fiR-lnJ In(l-j) ) + O (e') 
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+ 



u 



45 9 



2 - y - ^InnR - -ln^nj? + 41nj ln(l - j) 



(3.49) 



where 



nR^nRj{l-j) , (3.50) 

as above. The absence of any e-poles in Eq. (3.49) verifies the renormahzabihty of the 
persistence length to second order in u. 

To check the power law (3.35) we put ur — I, and we evaluate 

,1-21/ Lr (J, 1) 



(3.51) 



for u = u* in strict e-expansion, using Eqs. (3.12), (3.17). Some algebra yields 



y + f + -^(J)+-^(l-J) 



+ O (e^) 



(3.52) 



All the leading logarithmic singularities are eliminated, the remaining singular terms con- 
tained in (j) and {1 — j) being of order J In J and (1 — j) In (1 — j). They thus only give 
rise to subleading corrections to the dominant behavior L{j,n) ~ j^^'^-, j <^ n. This proves 
the asymptotic power law to second order in e. 

For a quantitative evaluation or our result we consider the ratio (3.39): 

L (j, n) 



2dPal (n) 



The second order result for a% {h) reads 



he 



(3.53) 



a 



E,R 



2 V 2 

/5 _ 

4 18 12 



InhR [l -'-)-'- In' hR + 0{e')) 
+ la,-hnhR + hn'hR + 0{e)^+0 {u') , (3.54) 



where 



"0 



dt 



\nt 



(3.55) 
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With relation (3.20) among L [j, n) and [j, ur) we find 



We now use Eqs. (3.49), (3.54) together with 



4 V 32 
to find in strict e-expansion 

f 21 

p{j,n) = l + -/ + _,V + -/(l-/)(21nn«-lnJln(l-j)) 

+|^/'(^-3+^ao + ^(j)+^(l-j))+0(6^) . (3.57) 
We first consider the fixed point / = 1, where p (j, n) reduces to 

P*(j) = l + | + ^(^ + Y + ^«0 + ^(j)+^(l-j))+O(e^) . (3.58) 

We note that the dependence on Inriji and thus on the choice of the renormahzed length 
scale has dropped out, as expected for a critical ratio at the excluded volume fixed point. 
To check the result (3.58), we may exploit the sum rule (1.3): 

n 

By virtue of the power law E?* {h) = (j(l — j))'^" R^* (n) , which holds in the excluded 
volume limit, the sum rule immediately takes the form 



1 



dj{j{l-j)f''-'p*{j)^l . (3.59) 



Employing the e-expansion (3.17) of u and writing 



P*{j)^l + l + e'p2{j) + 0{e') , (3.60) 



we find from Eq. (3.59) 



/ 



1 

2 



rfjP2(j) = ^ + ^« 0.1429 . (3.61) 





We have checked by numerical integration that the result (3.58) obeys this relation. 
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FIG. 3: Pq (j) (Eq. (3.64)) as function of j. Results of the e~expansion: O (e°), short dashes; 0{e), 
long dashes; O (e^), full line. One loop crossover model: dot-dashed line. Sum rule: dotted line. 

Despite its fairly complicated analytical form, (c.f. Eq. (3.48)), the function p2 (j) is 
almost independent of J. It monotonically decreases from p2 (0) = P2 (1) ~ 0.1496 to 
P2(l/2) ~ 0.1405. As a result, to the order considered p* (j) essentially can be taken as 
constant: 

p* (J) = 1.25 + p2(j)^p* = 1.393, (e = l) . (3.62) 

If p* (j) were strictly independent of Jto all orders of e, the sum rule (3.59) would yield the 
result 

P* = ^I^^ 1-408, (e = l) . (3.63) 

This is very close to the result of our second order calculation. The small difference, left for 
the higher orders, suggests that also in higher orders the J-dependence of p* (j) stays weak. 
To illustrate the J-dependence of L* (j, n) we in Fig. |21 have plotted the universal ratio 

Plij) = ^^^^^ = ijil-J)r-'P*{J) , (3.64) 

evaluated in the different approximations discussed in this section. Obviously the O (e^)- 
result (full line) is almost indistinguishable from the hypothetical result (dotted line) based 
on the sum rule p* (j) = p* ^ 1.408. The one loop crossover model (Eq. (3.31), dot-dashed 
line) coincides with the O (e^)-result in the tails, but deviates somewhat towards the center 
of the chain. This is understandable, since the choice of the reference chain length n and 
thus of the renormalized length scale is optimized for J or J ^ 1. However, also in the 
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center of the chain the one loop crossover model considerably improves the plain O (e)-result 
(long dashed line). 

So far we considered universal ratios. Turning to the normalized persistence length 
L [j, n), as defined in Eq. (1.1), we find from Eq. (3.64) 

L* (j,n) ^ , , R^* (n) 



For the discrete chain model bare perturbation theory shows that (s|) for d = 3 weakly 
depends on j. This is an endeffect, saturating for j ^ 1. Furthermore, taking the continuous 
chain limit we find 



2dil 



/,(/?e), {d>2) 



where the function fg is independent of J, provided < ] < 1. This suggests to introduce 
the arclength of the chain 



L,= ^{s])n , (3.66) 

neglecting any j-dependence. Employing the result (3.65) with p* (j) taken from the sum 
rule (Eq. 3.63), we thus find the simple expression 

L*{j,n)^l.m{jil-j)f^~'S_^^ (rf = 3) , (3.67) 

which generalizes a standard result |])| based on the wormlike chain model. It, however, must 
be noted that in the present context (s|) and thus Lc are effective nonuniversal quantities, 
so that Lc may differ from the arclength calculated from a physically realistic model by a 
factor of order 1. 

We finally consider the crossover from 0-conditions to the excluded volume limit. We 
use the representation 

p{j,n)aE{h) , (3.68) 



where a|; {n) is taken from Eq. (3.53): 



o^E [n) = —g- Oil^R {ur) 



Eqs. (3.14), (3.15), (3.29) yield 

\l-f\H^-^)H{f) . (3.69) 
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FIG. 4: L /6£^ as function of z = yj Snfi. Full line: O (e^); dashes: O (e); dot-dashed line: 

one loop crossover model. 

With the simple choice raij = 1 we find in strict e-expansion of a\ ^: 



ar, [n 



1 - 



Hif) 



(3.70) 



In Fig.lHwe plot the result for L (j = |, n) /2(i£^ as function of i = a/ s„ (/3e) a variable 
which is the counterpart in the renormalized theory of the standard 2;-variable of two- 
parameter theory It is related to the intermediate variable / via Eq. (3.15): 



l = /-2|l_/|l/-^^52 



(3.71) 



m if) 

For a comparison we also show the result of the one-loop crossover model (Eqs. (3.38), (3.37), 
u* = 0.364, no = 0.53). (We recall, that the choice ur = no = 0.53 is optimized to the one 
loop crossover model. No such analysis is available for the O (e^)-calculation, so that we 
used n/j = 1 as simplest possible choice.) Fig. IH illustrates the gradual increase of L{j,n) 
with increasing excluded volume strength. It furthermore again shows that the one-loop 
crossover model considerably improves the plain O (e)-result. 



IV. SIMULATIONS 



We measured the persistence length in the Domb- Joyce model, where the chain config- 
uration is modelled as a random walk on a regular lattice. Each configuration is weighted 
by a factor (1 — w)"^ /Z, where n2 is the number of pairwise intersections, and Z = Z (n) is 
the partition function. Using cubic or square lattices, we can identify a microscopic segment 
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FIG. 5: Simulation results for po{j,n) = nL{j,n) /B^{n) as function of J. Triangles: n = 20; 
crosses: n = 200; points: n = 2000. 



with a primitive lattice vector, which defines the unit of length so that L (j, n) = L (j, n). 
To simulate the model we used the PERM-algorithm developed by Grassberger in the 
form also employed in our previous work ^ on the correlations among segment directions. 



A. Results for d = 3 



On the cubic lattice our simulations extend to a maximal chain length rimax = 2000. Most 
measurements were performed for w = 0.4, which for this lattice is known to be close to the 
excluded volume value w*, where the leading corrections to scaling vanish, (0.4 < w* < 0.5, 
according to Ref. In renormalized variables w* corresponds to / = ^ = 1. To study 

the approach to the excluded volume limit we also used w = 0.01 (weak coupling: / < 1) 
and w = 1 (self-avoiding walks, strong coupling: / > 1). 

Results for the ratio po (j;"^) (Eq- (3.64)) with w = 0.4 are shown in Fig. |S| We note 
that the results are essentially independent of n, as expected for an universal ratio at the 
fixed point. Only for the shortest chain {n = 20) some small deviation from scaling can be 
identified in this plot. These indicate nonuniversal 1/ra-corrections. From the experimental 
side these results establish the scaling law. 

In Fig. |31 we have not included theoretical curves, since on the scale of that figure the 
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J 

FIG. 6: po (j) n) = fi L (j, n) /i?g (n) as function of j. Data: n = 200, large points; n = 2000, small 
points. Theory evaluated at the fixed point; full line: O (e^); dot-dashed line: one loop crossover 
model; dotted line: sum rule (3.63). 

O (e^)-result would just be covered by the data. For a comparison among theory and 
data we therefore magnified the plot by dividing out the power law j(l — j)^*^^^, result- 
ing in the ratio p* (j), (Eq. 3.39). To reduce the data we used the relation Rl{n) /h = 
(j(l — j))^*^^^ -Rg (n) /n, vahd in the excluded volume limit. Fig. IHl shows the results for 
n = 200 and n = 2000. We note some small effect of chain length n: the data for n = 2000 
seem to trace out a flatter curve, being lifted in the center and lowered in the wings com- 
pared to the data for n = 200. This is consistent with w = 0.4 being slightly below the fixed 
point value w*. For 0.2 < j < 0.8 the data essentially fall between the O (e^) result and the 
prediction of the sum rule. The increasing deviation among theory and data outside that 
range just indicates that reducing the data by dividing out the power law — j))^'^^^ 
strictly is adequate only at the fixed point and only for j ^ 1 and n — j ^ 1. With this 
in mind, and taking into account the scale of the figure, we may state excellent agreement 
among theory and data. For completeness we note that the scatter in the data for n = 2000 
decreases with increasing j. This is a feature of the PERM-algorithm, which within a single 
numerical experiment executes many correlated measurements. For the present problem the 
strength of the correlation depends on j. A fair impression of the statistical scatter is given 
by the data for small j. 

We finally consider the crossover towards excluded volume conditions. Fig. U\ shows 
a doubly logarithmic plot of L{n/2,n) /Qp as function oi z = ^/s^. The nonuniversal 
parameters i and v = y/s^ for the Monte Carlo model as used here have been determined 
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FIG. 7: logio /6^^ as function of lo 

Sio (^'^ ) = logio measured for w = 0.1, 0.4, 
1.0. Dot-dashed lines: one loop crossover model. The bars indicate nmax = 2000 for the different 
tt)-values. 

in previous work [3| by fitting quantities like the end-to-end distance to the one-loop 
crossover model. For consistency we therefore also here compare the data to the one-loop 
model (Eqs. (3.28), (3.71)), noting that the difference to the full O (e^) result is quite small, 
(see Fig. 1^. The theoretical curves show the well known two branched structure 0,^^ of 
the crossover scaling functions. The upper branch represents crossover towards 0-conditions 
(/ < 1), the lower branch corresponds to / > 1. This structure is nicely confirmed by the 
data. As mentioned above, in reducing the data we used the same parameter values £ and v 
as in Ref. and we allowed for an additional nonuniversal factor ~ 4/(8^), as discussed 



in the context of Eq. (3.67): independent of w we multiplied the data by 1.26. We should 
note that in principle we also should allow for corrections due to the discrete microstructure 
of the Monte Carlo chain. This is particularly relevant for weak excluded volume {w = 0.1), 
where the measured persistence length even for n = 2000 exceeds the microscopic bond 
length only by a factor of about 2. Also the prefactor ~ (s|) in principle should depend 
on w. Playing with such corrections we obviously can bring the data in the range of small 
z even closer to the theory. We, however, have not pursued this further, being content with 
the good overall agreement among theory and data exhibited in Fig. [7| 
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FIG. 8: p{j,n) as function of j, measured on the square lattice. Data: crosses, n = 200; points, 
n = 2000. Excluded volume parameter: w = 0.4. The curve gives a fit to Eq. (4.2). 



B. Some consideration of d = 2 

Grassberger and later Redner and PrivmanP, Q|, considered the persistence length 
on two dimensional lattices, using exact enumeration and simulations. They only discussed 
L (1, n), projecting the end-to-end vector on the direction of the first segment. They found 
that L (1, n) slowly increases with n, pointingto logarithmic behavior p: L (1, n) ~ Inn or 
to a power law with a very small exponent p|: L (l,r;,) ~ ^o.063_ ^^^y. ^^^^^ ^^^is is quite 
different from the behavior in three dimensions, where L {l,n) stays microscopically small: 
L {l,n) ~ a/ (sf), both according to our theory and to our and previous [71 simulations. 

For a closer examination of this effect we carried through simulations on the square lattice, 
using w = 0.4 and again going up to ^max = 2000. Results for the ratio 



N /-/I -^^ -1/2 nL{j,n) 



n] 



(4.1) 



are shown in Fig. |H1 (Note that 2z/ — 1 = 1/2 in = 2.) We again note a reasonable collapse 
of the data on a single master curve. Note that the fixed point value w* is unknown for 
d = 2, and w = 0.4 need not be close to w*. Thus the data collapse indicates that the chain 
lengths n > 200 in d = 2 are sufficient to reach the excluded volume limit. We thus believe 
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that our results indicate that the scahng law (3.34): 

also holds in two dimensions. However, the scaling function evidently does not obey the 
simple law 



-\\2u-l 



The difference to three dimensions is drastically illustrated by comparing Figs. 6 and 8, 
lich both up to a constant give C (j) 
are reasonably well fitted by the form 



which both up to a constant give C (j) / (j (1 — j))'^'^ ^. Indeed, in two dimensions the data 



P (j, n) 



1-b 



In(l-j) , In J 



(4.2) 



J 1-J 
which for J = 1/n yields logarithmic behavior 

L (1, 77.) ~ Inn , 

consistent with previous findings. This is illustrated by the full curve in Fig. |H1 which gives 
the expression (4.2) with b taken as fitparameter. a = a{b) was determined from the sum 
rule (3.59). The resulting parameter values are found as 6 = 0.275 and a = 1.376. 

Is there any support for the existence of such logarithmic terms from the side of the 
theory? It first must be noted that the e-expansion clearly does not yield such terms, but 
the naive extrapolation from d = 4 down to = 2 is somewhat doubtful. Indeed, evaluating 
the unrenormalized first order contribution (3.2), (3.3) directly for d = 2, we find a leading 
behavior 

/3A 0) « (i^^ + ^) . (4.3) 

(It is this result, which motivated the ansatz (4.2) for fitting C (j).) This logarithmic anomaly 
is specific to two dimensions and cannot be found in the e-expansion. In the context of 
polymer physics such anomalies first have been discussed by Des Cloizeaux. (See chapters 
12, Sect. 3.2.4 and 10, Sect. 4.2.6 of Ref. He found that such anomalies occur in 

dimensions d = 4 — 2/p {p integer). For the partition function Z (n) they yield a prefactor, 
which depends both on the usual variable z = Pen'^^'^ of the unrenormalized two parameter 
model and on In n. In standard observables like Rl this prefactor cancels and no anomalies 
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show up. The results of Refs. as well as those presented here suggest, that for 

the persistence length in d = 2 the anomaly survives in the form of terms depending on 
In (j/n) , In {{n — j) /n). This is not inconsistent with des Cloizeaux's findings, which only 
imply the absence of In ra-terms. Indeed, exploiting Eq. (12.3.49) of Ref. [l^ we immediately 
find that the reducible contribution L^^"^^^ (Eq. (2.9)) in bare perturbation theory for 

d = 2 picks up an anomalous prefactor 

'ln(l-j) , In J 



g{(3en,j) exp 



J 1 - J. 

where the constant Ci and the amplitude function g can be calculated perturbatively. Note 
that the combination of the anomalous terms is the same as in Eqs. (4.2) and (4.3), our 
fit formula (4.2) just taking into account the lowest order correction following from des 
Cloizeaux's work. 

For the irreducible contribution the same anomaly occurs. This is a necessary prerequi- 
site for renormalizability, which mixes L^^"^^^ and L'^"^^ (j, '^)- We, however, have not 
pursued the matter further, since a short calculation shows that also the normalizing factor 



(s^) of L (j, n) (Eq. (1.1)) is plagued by anomalies: (sp = 4£q(H- const [3^ In j) for n — oo, 
d = 2. This sheds some doubt on the very applicability of the self-repelling Gaussian chain 
model for a calculation of the persistence length in two dimensions. (We recall that (s^) in 
three dimensions rapidly tends to a constant.) Irrespective of this concern this discussion 
shows that also theoretical arguments support the existence of a logarithmic anomaly in two 
dimensions. 



V. CONCLUSIONS 

As expected, we have found that the persistence length of an excluded volume chain 
is a critical quantity, which can be evaluated within the framework of renormalized two 
parameter theory. We have proven the renormalizability to two loop order. Furthermore, 
considering the excluded volume limit of a long self repelling chain we both in theory and in 
simulations have found that the persistence length in three dimensions is well approximated 
by the surprisingly simple expression 

1.408 (l(l-i))"™S>), 
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where Lc is the effective contour length of the chain. This expression contains both the 
overall scaling L* ~ n^'^~^ and the dominant variation along the chain. The full O (e^)- 
result predicts some additional complicated dependence on J = j/n, but numerically this is 
completely irrelevant. Our result shows that in d = 3 the persistence length is of microscopic 
size for j close to a chain end, consistent with previous exact enumerations |7J. It rapidly 
increases towards the center of the chain. This is consistent with our previous results 8| 
on the correlations among segment directions, which show a similar end-effect. It may be 
of interest to give an example of the size of L* {j,n) for a typical system. For j = n/2 
and n = 2000 our simulations yield L* (n/2,n) ~ 3.36, to be compared to Re* (n) ^ 78, 
both measured in units of the lattice spacing, which here is identical to Kuhn's effective 
segment length. Since the chain length in the Monte Carlo model roughly corresponds to 
the polymerization index of a highly flexible polymer, this gives an impression of the typical 
persistence length in a medium size polymer coil. 

Outside the excluded volume limit our results show that L (j, n) essentially varies like the 
end-to-end swelling factor of a chain of effective length h = j {1 — j/n), with a prefactor 
which slowly increases from 1 (B-point) to about 1.4 (excluded volume limit). Again this 
result is in good accord with our simulations. 

For two-dimensional systems both our simulations and theoretical arguments point to 
the existence of a logarithmic anomaly: L{j,n) ~ Inn {n —>■ oo, j fixed), thus supporting 
previous findings . However, a precise analytical analysis might need a model with 

segments of fixed length, not the Gaussian chain model employed here. 

With respect to the much more delicate P| problem of an effective, locally defined persis- 
tence length in polyelectrolytes our results might be used to eliminate the excluded volume 
effects, if the persistence length is determined by fitting the observed coil radius to a worm- 
like chain model. This is a strategy sometimes followed in the analysis of data, (see e.g. 
Refs. 13, 13] )• Our results also point to a strong variation of the persistence length along 
the chain. 
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